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A microscopic theory for the coupling of intrinsic Josephson oscillations and dispersive phonon branches in 
layered superconductors is developed. Thereby the effect of phonons on the electronic c-axis transport enters 
through an effective longitudinal dielectric function. This coupling provides an explanation of recently observed 
subgap resonances in the Jdc-Vdc-curve of anisotropic cuprate superconductors forming a stack of short Josephson 
junctions. Due to the finite dispersion these resonances can appear at van-Hove-singularities of both optical and 
acoustical phonon branches, explaining low-voltage structures in the I-V-characteristic, which are not understood 
in phonon models without dispersion. In long junctions the dispersion of collective electron-phonon modes parallel 
to the layers is investigated. 



1. INTRODUCTION 

The main features of the electronic c- 
axis transport in the high-Tc-superconductors 
Tl2Ba2Ca2Cu30io+5 and Bi2Sr2CaCu208+5 be- 
low the critical temperature Tc can be well under- 
stood in the model of a stack of superconducting 
Cu02-layers, which form a onedimcnsional stack 
of intrinsic Josephson junctions [Q. 

One of the few phenomena, which are specific 
to the intrinsic JosephsonefFect and cannot be 
seen in arrays of conventional superconductors, 
is the interaction of the Josephson oscillations 
and phonons. This has been observed experi- 
mentally in the form of resonances in the I^c- 
Vdc-characteristics at low voltages ||^~|^]. Their 
main features could already be understood theo- 
retically in the framework of a simple lattice dy- 
namical model of damped local oscillators in the 
barrier 

Despite the clear evidence of the effect a gen- 
eral multi-band theory, in which parameters like 
the phonon frequencies or the oscillator strength 
in the dielectric function are specified microscopi- 
cally, is still missing and shall be presented in the 
following. This more general formalism is in prin- 
ciple appropriate to include the eigenvalues and 
eigenvectors of the dynamical matrix as obtained 



in any microscopic lattice dynamical calculation 
into the Josephson theory. Qualitative effects of 
the phonon dispersion are discussed in a simple 
rigid ion model, e.g. the van-Hove singularity at 
the upper edge of the acoustical band, which can 
be detected in the Jdc-Vdc-curve. Details of the 
derivation and further results can be found else- 
where 

2. THEORY OF PHONON COUPLING 

The Josephson effect between the layers n and 
n -f 1 is described by a RSJ-like equation 

i = ic.S^^ln{t)+ iqp{En{t)) + ^^K{t), 

where j is the bias current (density) , jc sin 7„ [t) 
the supercurrent, jqp{En) the quasiparticle cur- 
rent and "fn{i) the gauge-invariant phase differ- 
ence, which is related to the average electric field 
En{t) in c-direction in the barrier: h-jnit) = 
2edEn{t) {d: thickness of the barrier). In contrast 
to this, the displacement current density eo-E^(i) 
does not depend on En, but the field EP{t) set- 
up by the conduction-electron charge-fluctuations 
Spn on the layers, which obey the Poisson equa- 
tion 6pn{t) = eo{EP{t) - E^_^{t) and the con- 
tinuity equation j„(t) — jn-i{t) = —Spn{t). All 
quantities are constant along the layers (qn = 0) 
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as in short junctions. 

The crucial point is to determine the (hnear) 
relation between the field 

EP{q,,Oj) = E - E'°'' = €p^{q,,Oj)E{q,,Uj) (1) 

created by the charge fluctuations alone and the 
total electric field E{qz,uj) due to both electrons 
and ions. This can be described by the (longitu- 
dinal) phonon dielectric function 



where the suszeptibility (A: phonon branch) 



^ ui^iqzX) - 



(2) 



(3) 



and the oscillator strength 

|»feA)P^^z/-(-'^-")f;;'^-")z:,. (4) 

depend on the eigenvalues uj(qz, A) and eigenvec- 
tors e{K\qz\) of the dynamical matrix Dq,^(^^, ), 
which contains all the short and long range inter- 
action between the ionic cores in the insulating 
state. The eigenvalues u!{qz,X) therefore repre- 
sent the bare phononfrcquencies in the absence of 
the charge fluctuations dpn- The appearence of 
special q^-dependent effective charges Zf^{qz) in 
equ. ^ reflects the different contribution of ions 
on and between the superconducting layers and 
is one of the distinctions to the conventional di- 
electric function. The function epY^{qz,uj) has ze- 
ros at the longitudinal eigenfrequencies u;(qz,X) 
of modes with polarization in c-direction. 

A similar consideration 0| shows the general 
form 

x(g,w) 



eph(q,t^) = 1 + 



(5) 



1 - f{q}xiq,^ 

of the phononic dielectric constant, which reduces 
in special cases to the longitudinal Cp^iqz) — 
eph(9|| = 0, <7z) and transversal dielectric constant 



(9||) = eph{q\\,qz = 0) 



X{q\\,uj), which 



is relevant in optical experiments. 

With the knowledge of the phononic dielectric 
function eph the /dc-Vdc-curve can be calculated 
with the ansatz 



7nW 



'jQ + ujt + S^nit) n resistive 



70 + Sjnit) 



n else 



(6) 
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Figure 1. /dc-Vdc-curve of the first branch in the 
two- atomic chain model {Mi > Mb). 



by linearizing in the small oscillating parts 
Sjn{t) <C 1. The final result for the first branch 

(Fdc = hLu/{2e), Lol := 2edjc/{heQ)) 



jiV) 



■ ftn jc^p^ 1 



(7) 



exhibits a phonon resonance exactly at the zeros 
of the modified dielectric function 



1 



1 



92 ""ph 



-1 uj-" 



i.e. in the unrenormalized phonon bands uj{q, A). 
More specifically peaks develop at the van-Hove 
Singularities of the phonon density of states due 
to the averaging over qz in equ. ^. 

As a reliable microscopic lattice dynamical the- 
ory for the dynamical matrix is not yet avail- 
able and as the intrinsic Josepshon effect has to 
be treated in a phenomenological (one-band) pic- 
ture until the details of the pairing mechanism are 
understood, the above results will be illustrated 
qualitatively in a simple model of a two-atomic 
chain of masses M„ and charges Zi = —Z^. 

If the mass on the layer (Mi ) is larger than the 
one in the barrier (Afb), the acoustical van- Hove 
singularity is suppressed, but the optical band 
shows a double peak at the band edges (cf. fig. 
|l|), the distance being the bandwidth. The (more 
realistic) case of Mi < Mb exhibits single peaks 
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at the upper edges of both acoustical and optical 
bands and is compared with experiment in fig. |[ 
If several junctions are in the resistive state, 
their dynamics is coherent due to the coupling 
with phonons, even if no other (e.g. inductive) 
interaction is present. This is important for the 
application of short junctions as effective high- 
frequency devices. Analytical expressions similar 
to equ. for the second branch can in principle 
distinguish different phase locked solutions in the 
^dc-Kic-characteristic 

3. COLLECTIVE MODES 

The zeros of the real part of the total dielectric 
function 



etot(g,w) = eph(q,w) 



^ + — , (9) 



define the frequencies a;oi-ph(<f, 'i') of the collec- 
tive modes of the coupled conduction electrons 
and phonons, where the Josephson plasma fre- 
quency is given by {\ab'- coherence length): 



\2 

2-^{cos{q^d) 



1)' 



-)• (10) 



These modes are different from the bare phonon 
bands uj{q, A) and can in principle be observed in 
neutron scattering. In fig. |^ the plasma band 
u!pi{q} for ^ mixes with the longitudinal 
acoustical (LA) and optical (LO) phonon modes 
uj{q,X), which are almost dispersionless on this 
scale. Note especially that a frequency gap of 
width yj/2LUpiD, appears in the spectrum, which is 
formally (not physically) similar to the polariton 
dispersion of coupled phonons and photons 0. 

4. EXPERIMENTAL RESULTS 

Recently the explanation of the subgap reso- 
nances in Refs.[||-^ with the phonon coupling 
mechanism presented here could be well con- 
firmed by Raman measurements on the same 
samples [||J^ and infrared refiectivity experiments 
with grazing incidence (see table |^) . Note 

that in our theory also Raman-active or silent 
modes may couple to intrinsic Josephson oscil- 
lations for Qz ^ and that in contrast to optical 
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Figure 2. Dispersion cjoi-ph(9a;, <lz — 0.055) of col- 
lective plasma-phonon modes for <C Qz (hor- 
izontal lines: bare phonon frequencies uj{q,X), 

Aah/d= 100, Xj ^ dXc/2Xab). 



experiments (at the F-point) the whole BriouUin 
zone contributes in an average way (cf. equ. ||). 
This provides a reason for slight discrepancies of 
these experimental data in table ^ For more ex- 
perimental and theoretical references see 

The qualitative features of the subgap reso- 
nances can already be understood in the local 
oscillator model of ||^: The position of the reso- 
nance is completely independent on temperature, 
magnetic field or the geometry of the probe, while 



the intensity of the structure varies 



Jc 



with 



the critical current density jc(T, B). Also the be- 
haviour in external pressure is consistent with the 
phonon interpretation ||^. 

Beyond this, in the more general model de- 
veloped above resonances at van-Hove singular- 
ities, e.g. at the upper band edge of the acousti- 
cal phonon band, can be described. This might 
be an explanation for a peak seen in Q at 3.2 
mV (=1.54 THz) in the /dc-Vdc-characteristic of 
Tl2Ba2Ca2Cu30io, as the same frequency is pre- 
dicted in lattice dynamical calculations ||ll[ for 
the upper edge of the acoustical band and because 
there are no optical phonon bands at that low fre- 
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Table 1 

Comparison of the frequencies fsg — h Vdc/2e 
(in THz) of the most pronounced subgap res- 
onances and of infrared- {Jlo) f^nd Raman 
active {fro) modes in Bi2Sr2CaCu208 and 
Tl2Ba2Ca2Cu„02„+4- 



Phonons in Bi2Sr2CaCu208 


fsg 2.97 3.89 5.17 5.60 | 


5 


1 


/lo 2.85 5.07 1 


9 


1 


/lo 2.86 5.16 




m 


/to 3.80 






= 1 


Phonons in Tl2Ba2Ca2Cu„02n 


+4 






/,g 3.63 4.64 1 


5 


n 


= 3 


/lo 4.50 j 


9 


\n 


= 2 



quency (cf. Fig. ^). Also the double peaks shown 
in fig. 1^ might have been seen in the satellite 
structures at 5.17 THz (10.7 mV) and 5.6 (11.6 
mV) in the /dc-Vdc-curve of Bi2Sr2CaCu208+5 
Q , which is consistent with the theoretical width 
- 0.3 THz of this phonon band Q. 

5. CONCLUSIONS 

In this contribution the microscopic theory for 
the coupling between Josephson oscillations and 
longitudinal phonons in intrinsic Josephson sys- 
tems like the highly anisotropic cuprate super- 
conductors has been developed. The analytical 
result equ. (^ for the /dc-Vdc-curve for one re- 
sistive junction and the detailed form of the lon- 
gitudinal dielectric function equ. (||) describing 
the coupling has been obtained. In paritcular, 
not only optical but also acoustical phonons at 
the edge of the Brillouin zone couple to Joseph- 
son oscillations explaining a structure observed 
in 0] in the /dc-Vdc-curve occurring at low volt- 
age. The discussion of the dispersion of collective 
phonon-electron modes is an important first step 
to understand the effect of phonons in long junc- 
tions 0. 
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W-7405-ENG-36 (C.H.) for financial support. 
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Figure 3. Fit of the experimental /dc-Vdc-curve 
in Tl2Ba2Ca2Cu30io Q near the band edge of 
the acoustical branch (at 1.5 THz) and an optical 
branch within the two-atomic chain model (inset: 
first two resistive branches). 
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